Diffuse first-order piiase transition in NaNbOsiGd 
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Although diffuse phase transitions have been Icnown 
for a long time, in ferroelectric relaxors 0, 01, they 
have not been studied as much in perovskite-based an- 
tiferroelectrics |^ ^ Despite this, the latter have 
some obvious advantages in technical applications be- 
cause of the absence of toxic lead and strong frequency 
dispersion of dielectric permittivity, e', inherent to relax- 
ors. NaNb03:Gd(NNG) crystals were grown by the flux 
method At a; < 0.1 e' has hysteresis (hereafter 
we mean thermal hysteresis of e'), which disappears at 
higher x (Fig. 1). Fig. 2 (a-b) shows e' obtained at 
heating of a NNG9 crystal up to T = Ta with subse- 
quent cooling (Fig. 2a), and cooling the crystal down to 
Tcooi with subsequent heating (Fig 2b). It is seen that 
the hysteresis loop area depends on Ta and Tcooi strongly 
that implies that the diffuse first order phase transition 
in NNG9 develops step by step in finite volumes. As a 
result, the crystal's state becomes inhomogeneous. The 
problem of describing dielectric permittivity in inhomo- 
geneous media is usually discussed in terms of the effec- 
tive medium approximation 9] , within which we obtained 
that the fraction of the low-temperature (LT) phase is: 



n(T) 



[e,{T)-e'{T)][e,{T) + 2e'{T)] 
3e'(T)[e2(r)-£i(T)] 



(1) 
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FIG. 1: Dielectric permittivity in NNGO (1), NNG9 (2) and 
NNG12 (3) measured at 100 KHz. The number after NNG 
shows the content of Gdi/aNbOa in mole percent. 
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FIG. 2: Experimental data on dielectric permittivity of 
NNG9 obtained (a) for different annealing temperatures: 
Ta \°C] = 150 (1), 170 (2), 190 (3); (b) different cooling tem- 
peratures: Tcooi [°C] = 90 (4), 110 (5), 130 (6), 150 (7), 180 
(8); c) The LT-phase fraction, n, obtained from experiment 
(squares); Solid lines are guides for the eye; (d) Comparison 
of the fit of expression ^ (solid line) to experimental data 
(circles). Dotted lines show normal first-order phase tran- 
sition behaviour. The angles between the dashed and solid 
lines characterize the degree of the diffuseness of the phase 
transition. 



where £i{T) and £2{T) are the dielectric permittivi- 
ties in the LT- and high-temperature (HT) phases, re- 
spectively. All these quantities can be obtained from the 
experiment performed. In agreement with our initial as- 
sumption, the n{T) dependence (Fig. 2c) is diffuse. One 
can also use the hysteresis loop area as a (nonlinear) mea- 
sure of n{T) on cooling and 1 — n{T) on heating. Below, 
we will discuss the shape of n{T) in more details. 

Temperature-dependent optical studies of NNG have 
been carried out by the rotating polarizer method, us- 
ing the Metripol (www.metripol.com) microscope sys- 
tem and a precise heating stage (Linkam HFS91). 
NNG crystals with x = 0.02 - 0.09 display distinctive 
changes of birefringence with temperature similar to that 
in NaNbOs. At the phase transition point, a spontaneous 
splitting into small regions (less than 0.001 mm in size) 
with diffuse boundaries occurs, and the distribution of 
the birefringence image becomes very complex. Afar of 
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FIG. 3: Theoretical curves for the model free energy (a): 1. 
(2.4, 0.1); 2. (2.4, 1.0); 3. (2,4, 1.5); and order parameter 
n (b-c): 4. (2.4, 0.1); 5. (2.4, 1.0); 6. (3.0, 0.1); 7. (1.5, 
0.1) (The values in parentheses are 'yz/ksT and a/kBT). 

the phase transition temperature, the scale of this non- 
uniformity increases up to 0.05 mm. 

The explanation of the experimental results obtained 
can be given in the same way as for relaxors P, . One 
can introduce a distribution function for the Curie tem- 
peratures, /(Tc), which, in the simplest case, can be 
described by a Gaussian function (this distribution is 
caused by internal local fields and stresses introduced by 
impurities 

s{T) = eo{T) + biT) j f{T,)e{T -T,)dT,^ 

- eo + Serf[(r-Teo)/r] (2) 

where £q{T) and B{T) are monotonic functions of tem- 
perature, which, in the first approximation, can be given 
by linear functions. Note that the Curie temperature 
of the ferroelectric phase transition is sufficiently lower 
than the temperature of the step. 9{x) is a step func- 
tion, which equals at x < and 1 at a: > 0; erf(x) is the 
error function. The fit (Fig. 2d) shows that the width 
of the distribution function for NNG9 is about 27 K on 
heating and 35 K on cooling, which is nearly comparable 
with the hysteresis width, 44 K. We found that the dis- 
tribution function width decreases with decrease of the 
Gd content. The results obtained allow one to suppose 
that the disappearance of the hysteresis and the dramatic 
increase of the e(T) diffuseness observed experimentally 
in NNG12 (see Fig. 1), may be due to a crossover 
between the widths of the £(T) hysteresis and the Curie 
temperature distribution function. 

The diffuseness of a first order phase transition means 
coexisting two phases. This must increase the surface 
tension energy. Below, we discuss this surface tension 



contribution within a microscopic approach, which is an 
extension of a lattice gas model [loj : 

H — A ^ Hi + 7 ^ [ui + Uj - 2n.inj] (3) 

Here A is a temperature dependent difference in chem- 
ical potentials of the phases, A — k {T — Tc), where we 
suggest Tc to be given randomly with a Gaussian distri- 
bution function at T = TcoJ the eigen values of rii are 
(HT phase) or 1 (LT phase); in the the second term we 
introduce the interface energy due to surface tension (it 
is positive when the nearest sites are different and zero 
when they coincide) . The summation runs over the sites 
each of which is occupied by a solid either in the LT or 
HT phases. The summation in the last term is over the 
nearest neighbors only. We have derived the free energy 
(at a fixed value of A), in the mean field approxima- 
tion: F(A) -/fcsTln (l + e^/^"'^) + jzn'^ where h = 
—A — 2z7(0.5 — n) is a mean field and n —< Ui > is the 
average occupation probability of the LT phase. F and 
n may be averaged with a Gaussian function of A/ksT 
having the width a. We obtained that, the free energy 
has two minima at a small a and at 7 > IkoT j z where 
z is the number of the nearest neighbors (Fig. 3a). In- 
creasing the width of the distribution function decreases 
the barrier height between the states and decreases the 
inclination of the n(A) curve (the measure of the sharp- 
ness of the phase transition). The real situation is, of 
course, more complex because of the random distribu- 
tion of barriers in space (the decrease of the barriers may 
result in percolation of the new phase as it was found for 
manganates but the model does give a reasonable 

explanation of observed behaviour: the phase transition 
is diffuse in spite of the surface tension the role of which 
is suppressed due to spreading the Curie temperatures. 
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